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ABSTRACT 

The interplanetary magnetic fluctuation spectrum obeys a Kolmogorovian 
power law at scales above the proton inertial length and gyroradius which is well 
regarded as an inertial range. Below these scales a power law index around —2.5 is 
often measured and associated to nonlinear dispersive processes. Recent observa- 
tions reveal a third region at scales below the electron inertial length. This region 
is characterized by a steeper spectrum that some refer to it as the dissipation 
range. We investigate this range of scales in the electron magnetohydrodynamic 
approximation and derive an exact and universal law for a third-order structure 
function. This law can predict a magnetic fluctuation spectrum with an index of 
— 11/3 which is in agreement with the observed spectrum at the smallest scales. 
We conclude on the possible existence of a third turbulence regime in the solar 
wind instead of a dissipation range as recently postulated. 

Subject headings: magnetic field — MHD — solar wind — turbulence 



Introduction 



Turbulence plays a central role in a wide ra nge of astrophysi cal plasmas. Examples ar e 
given by the solar wind (IMatthaeus et al.lll999l). the inters tellar (jScalo &: Elmegreenl 12004 ). 
galactic and even intergalactic media (IGovoni et al.ll2006l ). In the solar wind, turbulence 
evolves freely and is not perturbed by in situ diagnostics, therefore it provides an ideal 
laboratory for studying high Reynolds number plasma turbulence. Thi s unique situation 
allows us to inves t igate for example the origin of aniso tropy (see e.g. iKlein et al.l Il993l : 
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( IKiyani et al.l 120091 ). or to analyze different regimes of turbulence characterized by a steep- 
ening of the magnetic field fluctuations spectrum with a power law index going from —5/3, 
at frequencies low er than 1Hz, to indices lying around —2.5 at higher frequencies (see e.g. 
Smith et alJbood ). 



The spectral break near 1Hz has been a subject for intensive studies and controversies 
in the last decades. It was first i nterpreted as the on set of dissipation caused for example 
by kinetic Alfven wave damping (ILeamon et al.lll998l ). Then, it was demonstrated that the 
wave damping rate usually increases very strongly with wavenumbers a nd should lead to a 
strong cutoff in the power spectra rather than a steepened power law (ILi et al.l 120011 ). In 
the meantime, there are some indications that the fluctuations are accompanied by a bias of 
the polarization sugge sting; the presence of rig ht-hand polarized, ou tward propagating waves 
( [Goldstein et al.lll994i ) . Also it was proposed (jStawicki et al.ll200ll ) that Alfven - left circu- 
larly polarized - fluctuations are suppressed by proton cyclotron damping and that the high 



frequ ency power law spectra are likely to consist of whistler fluctuations ( IMatthaeus et al. 
20081 ). It is currently believed that the steepening of the spectra at 1Hz is mainly due 



to non-linear dispersive p rocesses that range from kinetic Alfven waves (H asegawa fc Chen 



19761 ; iHowes et al.l 120081 ) , electr omagnetic ion-cycl o tron Alfven waves (IGarv et al 



or /and electron whis tler waves (IGhosh et al.l Il996t iGaltierl l2006t iGaltier fc Buchlin 



2008 



2007 



Narita &: Gary! I201CH ) in the framework of Hall magnetohydrodynamics (MHD) or simply 
electron MHD. 

The most recent solar wi nd observations made with the high resolutio n magnetic field 
data of the Cluster spacecraft (jSahraoui et al.ll2009t lAlexandrova et al.ll2009h reveal the pres- 
ence of a third region - called dissipation range - at scales smaller than d e and characterized 
by even steeper magnetic fluctuations spectra with a power law index arou nd —3.8. Thes e 



spect ra observed only on half a decade are interpreted as either a power law (jSahraoui et al. 



2009 ) or an exponential law ( Alexandrova et al. 2009h . Alth ough the theoretical interpre 



tation of such a regime is still open (IMatthaeus et al.l 120081 ) . a recent theoretical analysis 
shows that a kinetic Alfven wave casca de subject to collisio nless damping cannot reach elec- 
tron scales in the solar wind at 1 AU (IPodesta et al.l |2010| ) . The direct consequence is that 
the spectra observed must be supported by another type of wave modes. It is noteworthy 
that this new regime at electron scales gives rise to the same controversy as the steepening 
found two decades ago around 1Hz which brings up naturally the following question: Have 
we really found the dissipation scale of the solar wind plasma or is it the onset of a new 
turbulence regime? 



In this article, we investigate the turbulence regime at scales smaller than the electron 
inertial length d e through the electron MHD approximation. The assumption of homogeneity 
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and isotropy will be made to derive an exact and universal law for third-order structure 
functions. We show that this law corresponds to a magnetic fluctuations isotropic spectrum 
in /c -11 / 3 compatible with the solar wind measurements. Although the assumption of isotropy 
is in apparent contradiction with the observations, it is claimed that the method used is a 
powerful way to have a first estimate of the anisotropic spectrum. Indeed, the main source 
of anisotropy is the presence of a large scale magnetic field which reduces the nonlinear 
transfer along its direction. Then, the most relevant spectral scaling is the transverse one for 
which the spectral index corresponds to the isotropic case if arguments based on the critical 
balance condition are used. It is only in the asymptotic limit of wave turbulence - for which 
anisotropy is stron g - that the spe ctral index for transverse fluctuations is (slightly) modified 
(see the review by iGaltierl ( I2009a[ )). Finally, we conclude the paper on the possible existence 
of a third inertial range for solar wind turbulence instead of a dissipation range as recently 
postulated. 



2. Electron magnetohydrodynamics 



Electron MHD provides a fluid description of the plasma behavior on length scales 
smaller than the ion iner tial length (L and on time scales of the order of, or shorter than, 
the ion cyclotron period ( IKinsep et al.lll990f ). In this case ions do not have time to respond 
because of their heavy mass and merely provide a neutralizing background. Then, the 
plasma dynamics is governed by electron flows and their self-consistent magnetic field. This 
model has attracted a lot of interest because of its potential applications in fast switches, 
Z-pinche s, impulsive magnetospheric /solar corona reconnection and ionospheric phenomena 



^see e.g. 



Bhattacharied 12004 ; IChacon et all 120071 ). 



The inviscid three-dimensional electron MHD equations can be written in SI as (IBiskamp et al. 



1996|) 



d t (l - d 2 e A)B = -diV x [J x (1 - e^A)B] , (1) 

where B is the magnetic field normalized to a velocity (B — > JJiqP B) and J = V x B 
is the normalized current density. Under the limit of electron MHD we remind that the 
current density is pr oportional to the electron velocity. This equation has two invariants 
(IBiskamp et al.l 119991 ) which are the total energy 



E = \ J{B 2 + d 2 e J 2 )d*. 



and the generalized helicity 



H = / (A - d 2 e 3) ■ (B - <%AB)dx . 



(2) 



(3) 
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with A the normalized magnetic potential. 



Equation is often used when d e — > 0, namely for scales between di and d e , for 
which a /c~ 7 / 3 -isotrqpic magnetic energy spectru m is found numerically and heuristically 
(IBiskamp et al.lll996l Il999l ; iDastgeer et al.ll2000al ). This result i s compatible w ith a rigorous 
derivation of a universal law for third-order correlation tensors ( Galtier 2008al lbl). A steeper 
magnetic spec trum in fc~ n / 3 may also be found when the kinetic energy overtakes the mag- 
netic energy (IGaltier fc Buchlinl 120071 ) . Such a situation - generally not discussed in the 
literature - can only be observed when the full Hall MHD system is considered. Note that 
this —11/3 power law index, valid for length scales larger than d e , has a different origin from 
the one derived in the present paper which is applicable for scales shorter than d e . 



The behavior at scales s horter than d e has attracted much less attention (IBiskamp et al. 



19961 ; IDastgeer et al.ll2000bl ). This regime corresponds to the limit d e — > +oo for which we 
have <9fAB = — rfjV x [J x AB], or equivalently with the relation AB = —V x J, 



di 



-Jx(VxJ)-V$, 



(4) 



where $ is an unknown function. The second term in the right hand side may be seen as a 
gauge; actually, an analysis performed directly on the generalized Ohm's law shows that it 
corresponds to an electron pressure. Note that the form of equation (T4J) is well adapted to 
the problem under consideration since we are going to assume isotropy which means we will 
not consider any background (large scale) magnetic field B . 



3. Universal law for r < d e 

In the following, we shall derive an exact and universal law for third-order structure 
functions for homogeneous three-dimensional isotropic electron MHD turbulence at scales 
smaller than d e and discuss the implications in terms of magnetic fluctuations spectrum. 
After simple manipulations on equation (T4|) we obtain for the zth-component 

\d t J, = J e d e J t - d t ($ + J 2 /2) , (5) 

where the Einstein notation is used. Note that we also have dgJe = 0. The second term in 
the right hand side is similar to a pressure whereas the first term exhibits a sign of difference 
with the usual advection term encounters in Navier-Stokes equations. Actually, the Navier- 
Stokes equations may be recovered when the electron velocity is used instead of the current 
density or when the generalized Ohm's law is directly used. We made the choice to use the 
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well-known electron MHD equations (JT]) and (jSJ) mainly because we shall compare eventually 
the new prediction with the previous one for scales larger than d e (IGaltierl l2008aj ) . 



It is straightforward to derive a universal law for an homogeneous and isotropic electron 
MHD turbulence. First we introduce the second-order correlation tensor 



i2y(r) = (J,(x) J^x')) = (J t J') 



where x' = x + r. We obtain 



3 1 ' 



(6) 



—d t Rij{r) 



= {Jitytf) - + J* /2)) 

+ (J'JedeJi) - (J'di($ + J 2 12)) . 



(7) 



After simple manipulations where we use the homogeneity and the divergence free condition 
we get 



-rd t Rij 

di 



d n ((JM) - (J'M)) . 



Note that the pressure- type contributions removed as usual for isotropic turbulence (IBatchelor 
19531 ). When the diagonal part of the energy tensor is only retained we have 



d 



d t R 



t 1 Hi 



—2d rt (J i J^J i ) 
-2V • (3(J t Ji)) . 



(9) 
(10) 



At this level of analysis it is necessary to say a word about the small scale dissipation and 
large scale forcing terms which have been neglected so far. The dissipation is a linear term 
which is seen as a sink for the energy. Since we are interested in a universal behavior of 
turbulence we are in a situation where the scales considered are supposed to be much larger 
than the dissipation scales: in other words, we are deep inside the inertial range where the 
dissipation has no effect. The forcing term is assumed to be at the largest scales and acts 
as a constant source of energy for the system. Formally, the introduction of a small scale 
dissipation V and a large scale force T leads to the expression 



-d e d t Ru 



-dtd 2 e V 



(j(JiJl)) + r+v. 



An exact relation may be derived for third-order struc ture functions by assuming the 



following; assumptions spe cific to fully developed turbulence (jKolmogorovlll94ll ; lFrischlll995 



Politano &: Pouquetlll998l ). First, we take the long time limit for which a stationary state is 
reached with a finite mean energy dissipation rate per unit mass. Second, we consider the 
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infinite magnetic Reynolds number limit for whic h the mean energy di ssipation rate per unit 



mass tends to a finite positive limit, e J (see e.g. iBiskamp et al.lll996l ). By noting that 



Rn = (J 2 )-^(5J t 5J l ), (12) 

where 53 = J(x + r) — J(x), we obtain 

^dldtRa = d t {^d 2 e J 2 ) - ^dedtiSJiSJi) , (13) 

where the first term in the right hand side is the time variation of energy. Therefore, in 
the stationary state both terms in the right hand side are equal to zero. Since dissipation 
effects are negligible we only have to include the mean energy injection rate per unit mass 
e J . It is important to remind that the energy (an inviscid invariant) is built directly from 
the current density hence the name e J . As it will be shown below, this remark turns out to 
be fundamental for the prediction of the magnetic fluctuations spectrum. The insertion of 
the previous statements into (TTT1) leads to 



did 2 e V- (J(JiJl)) =e J . (14) 
The introduction of structure functions for the current density gives event uaH)!^ 

did 2 e V ■ {53(8 J) 2 ) = 4e J . (15) 



An integration of ( I15p over a full sphere of radius r (since isotropy is assumed) and 
the application of the divergence theorem give finally the universal and exact law for three- 
dimensional homogeneous isotropic electron MHD turbulence for scales smaller than d e ; it 
writes 

d l dl(5J L {53f) = -e J r 1 (16) 

where L means the longitudinal component of the vector, i. e. the one along the direction r. 
Note the positive sign in the right hand side which is compatible with the negative sign in 
front of the nonlinear term of equation (jlj). The most remarkable aspect of this law is that 
it not only provides a linear scaling for third-order structure functions within the inertial 
range of length scales, but it also fixes the value of the numerical factor appearing in front 
of the scaling relation. 



1 It is also possible to consider a system without external forcing (jLandau fc Lifchita ll989). In this case, 
we have to deal with the decay problem for which the time derivative of the energy is equal (up to a sign) 
to the mean energy dissipation rate per unit mass e J . It is still possible to assume the time independence of 
the second term in the right hand side of equation (1131) and to finally recover the same relation as (|15[) . 



-7- 



Extension to r > gL 



For scales larger than d e (but still smaller than dj) the universa l law for three-dimensional 
isotropic electron MHD turbulence takes the form (iGaltierl l2008al ) 



^([(JxB)xB'] L ) = -Vr. 



(17) 



where L still means the longitudinal component. This law may also be written as (IGaltier 
2009bh 

di<[(J3<B) x 6B} L ) = - l -e M r , (18) 

where X = (X + X')/2. Both universal laws f[T5]) and f[T8l may be gathered by noting that 
in this case the forcing scale is pushed at scales much larger than d e . One needs to consider 
the following expression in the stationary state and in the inertial range 

1 



-{d 2 e d t Ru + d t R ti ) = AfC + T + V = ATC + e 1 



(19) 



where Rjj = (BjB'j), MC is the nonlinear contribution and e T is the mean total energy in- 
jection rate per unit mass. Then, one obtains the general law for three-dimensional isotropic 
electron MHD turbulence (with r < dj) 

4 



4di([(J x B) x SB] L ) - djd 2 e ({53)' 2 5J L ) = ~^ T r. 



(20) 



This universal law conserves the linearity in r but emphasizes the role of each term according 
to the scale considered. The importance of each term is given by the current density which 
involves a derivative: at small scales (scales smaller than d e ) the contribution of the current 
density will be more pronounced than at large scales (scales greater than d e ). 



5. Magnetic fluctuations spectrum 

The universal law ffTB]) gives a precise description of electron MHD turbulence at scales 
smaller than d e . It also provides the possibility to predict the form of the magnetic fluc- 
tuations spectrum and then a comparison with observations. Dimensionally relation (fl6|) 
corresponds to 

did 2 e J 3 ~ £ J r , (21) 

which means an energy spectrum in 

/ F J \ 2/3 

Ej{k) ~{M) " ' (22) 
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comp atible with the dynamical equation (jSJ) and direct numerical simulations ( iBiskamp et al. 



19961 ). Since the current density and the magnetic field satisfy the equation J = V x B, we 



obtain the scaling relation 



n3 

d i( P—~£ J r. (23) 



1 e r 3 



Then, the corresponding magnetic fluctuations power spectrum scales as 

*<*>~G&r*~ u/ '- (24) 

This magnetic spectrum is significantly steeper than the energy spectrum. It is compatible 
with solar wind measurements where the magnetic field fluctuations are generally used to 
investigate turbulence at small scales, i. e. for freque ncies / higher than 1Hz. A power law 



around / has been reported (ISahraoui et al.ll2009l ) which is significantly steeper than what 



it is found at smaller frequencies where we have a spectrum around f~ 2 ' 5 . Note that the 
linear law in r ( 120]) corresponds in fact to a double-scaling law for the magnetic fluctuations 
spectrum in k~ 7 ^ 3 and k~ n ^ 3 for, respectively, length scales larger and shorter than d e . If we 
follow the Taylor frozen-in-flow hypothesis - which is questionable at this length scales - then 
we arrive at the conclusion that the general law (|20|) , although its relative simplicity with the 
assumptions of isotropy and incompressibility, is in good agreement with the observations. 



Conclusion 



The turbulence regime at scales smaller than the electron inertial length d e has been 
investigated through the approximation of electron MHD. A new universal and exact law has 
been established in terms of structure functions for the current density. This law leads to 
the prediction of a A; -11 / 3 power law spectrum for the magnetic field fluctuations compatible 
with the most recent observations made with Cluster. It is proposed that electron MHD 
turbulence provides a valuable first order approximation for the solar wind dynamics in 
particular below the length scale d e . It also provides the first prediction for the magnetic 
fluctuations spectrum at these length scales. The possibility to get a turbulence regime at 
electron scales questions the origin o f dissipation in the solar wind and more gene rally in 
space plasmas. In previous analyses ( ISahraoui et al.ll2009t lAlexandrova et al.l 120091 ) it was 
suggested that the range of scales where the heating occurs was discovered but it was also 
confessed that the characteristics of turbulence in the vicinity of the kinetic plasma scales 
are not well known neith er experimentally nor theoretically and are a matter of debate (see 
also iPodesta et al.l 120101 ) . It is believed that the present theoretical prediction will help 
significantly in such a debate. 
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A. Extension to 2D or slab turbulence 



We shall discuss now the extension of the exact prediction (1161) to the non isotropic case. 
The simplest situation is when a background magnetic field B is applied to the plasma flow 
for which the dynamic s becomes statistically axisymmetric. We will adopt the view point of 
MacBride et al.l ( 120081 ) and assume by simplicity that electron MHD turbulence is either 2D 
or slab. If we define the energy flux vector F as (see relation (fT5|) ) 



F = d i( i 2 e V ■ (6J(6J) 2 ) , 

then we obtain the general formulation for isotropic electron MHD turbulence 

4 



Fir) 



D 



e J r , 



(Al) 



(A2) 



where D is the space dimension. 2D turbulence means that F and B are perpendicular 
whereas they are parallel for slab (or ID) turbulence. Under these simplifications it is 
straightforward to derive the following predictions 



and 



F 2D (r) 



nsiafe / 



2£ Vswt sin 9 . 



4e Vswt cos 9 . 



(A3) 



(A4) 



where r is the time, and 9 is the angle between the mean magnetic field and the solar 
wind velocity V$w directions. Note t hat to obtain these predictions the Taylor hypothesis 
has been used ( MacBride et al.l 120081 ). It is believed that such relations could be useful for 
analyzing solar wind turbulence and to evaluate the perpendicular and parallel mean energy 
dissipation rates per unit mass. Note that a more sophisticated approach has been recently 
used to describ e axisymmetric electron MHD turbulence for length scales larger than d e 
( lGaltierll2009bl ). Under the assumption of critical balance it leads to a vectorial relation for 
the energy flux with a dependence in both r and 9. 
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